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Berry curvature and 4-dimensional monopole in relativistic chiral kinetic equation 
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We derive a relativistic chiral kinetic equation with manifest Lorentz covariance and a 4- 
dimensional Euclidean Berry monopole. The theory is based on Wigner functions of spin- 1/2 mass- 
less fermions in a constant electromagnetic background. By integrating out the po-component of 
the 4-momentum p, the previous 3-dimensional results derived (without vorticity) from the Hamil- 
tonian approach is reproduced. The phase space continuity equation has a source term proportional 
to the product of the vorticity and electric field while the axial anomaly arises from the flux of the 
monopole. This makes the chiral magnetic effect, vorticity, chiral anomaly, Berry phase, and the 
monopole can all be described in a unified way by Wigner functions. 

PACS numbers: 25.75.Nq, 12.38.Mh, 13.88.+e 



Introduction. — Chiral anomaly is an important quan- 
tum effect which is absent at the classical level. Recently 
it has been shown that such a microscopic quantum ef- 
fect can have a macroscopic impact on the dynamics of 
relativistic fluids. The chiral magnetic and vortical ef- 
fect (CME and CVE) are manifested as currents 
induced by magnetic field and vorticity. Such effects 
and related topics have been investigated within a va- 
riety of approaches, such as AdS/CFT correspondence 
[IHij], relativistic hydrodynamics j^-Hlj. and quantum 
field theory @,[ll-[l7|. Derivation of a semi-classical chi- 
ral kinetic equation incorporating CME and CVE will 
deepen our understanding of the underlying mechanism 
of these effects. Recently it has been found that fea- 
tures of Berry phase due to a 3-dimensional momentum 
monopole emerge in a chiral kinetic equation without 
manifest Lorentz covariance 



and 4-dimensional Euclidean monopole. These results 
reveal the inherent connection between the Berry phase 
and gauge invariant Wigner functions. One advantage of 
our approach is that the vorticity effect can be derived 
straightforwardly which, apparently, is not the case in 
other approaches. The non-covariant version of our co- 
variant equation is shown to be in agreement with previ- 
ous results if vorticity terms are turned off 3, 19| . Our 
derivation is qui te gene ral and valid not only for Fermi 
liquid as in Ref. 



is, m 



A semi-classical 
kinetic equation has also been derived in an electron sys- 
tem with Berry curvature [2fjj |. The Berry phase is a 
topological phase acquired by an energy eigenstate af- 
ter it undergoes an adiabatic evolution along a loop in 
parameter space 2l|. The Aharonov-Bohm effect is an 
example of the Berry phase in clcctromagnetism. The 
Berry phase has been applied to a variety of condensed 
matter phenomena, such as the quantum Hall effect, the 
spin Hall effect, etc. For a recent review, see e.g. Ref. 

In this paper we will derive a new chiral kinetic equa- 
tion with manifest Lorentz covariance from the Wigner 
function [2J|. We will show that such a chiral kinetic 
equation incorporates features of the Berry curvature 



19J but for any relativistic fermionic 
system. The phase space continuity can be shown to be 
broken by an anomalous term proportional to the prod- 
uct of vorticity and electric field. So the phase space 
measure is not conserved. It is modified by a factor re- 
lated to the Berry curvature. We will also show that 
the conservation laws of the right- and left-hand currents 
is broken by anomalous terms, which can be given by 
the flux of a 4-dimensional monopole in Euclidean mo- 
mentum space. Therefore a variety of properties such 
as CME/ CVE, chiral anomaly, Berry curvature and 4-d 
Euclidean monopole can be described in a unified way in 
the quantum kinetic theory with Wigner functions. 

Equation of motion with Berry curvature in 3- 
dimensions. — We will follow an example of Ref. [l8} to 
illustrate the concept of Berry curvature. We consider a 
Hamiltonian H ' = a ■ p for spin- 1/2 fermions in addition 
to the normal part H(p, x), where <x are Pauli matrices. 
Under an adiabatic evolution, the action of the helicity 
plus state in path integral is 

S= ( dt(x.p + x.A(x)-p-a(p)-H(p,x)], (1) 
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where A(x) is the electromagnetic vector potential and 
a(p) is the vector potential in momentum space resulting 
from diagonalizing H' in path integral. We can generalize 
the coordinate variables by combining p and x, £ a = 
(p, x) with a = 1,2, •••,6. The action can be cast into a 
compact form, 



s = J dt[- 7a (oL - m% 



(2) 



quantum interactions and we will show that all currents 
including chiral anomaly can be derived from the above 
equation. The Wigner function can be decomposed in 
terms of 16 independent generators of the Clifford alge- 
bra whose coefficients are scalar, pseudo-scalar, vector, 
axial vector and tensor respectively. The quantum ki- 
netic equation for the Wigner function then leads to two 
decoupled sets of equations for its components [2514271 ] . 
one of which relevant to our study reads, 



where 7 a = [a(p), — p — A(x)]. The equation of motions 
follows the Lagrangian equations, 



(3) 



0, p^ = 0, (6) 
0, VW M = 0, (7) 

(8) 

-2^^-p^), (9) 



where j ab = <9 a 7b(£) ~ <Vy<i(0 and given by 
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(4) 



where fl = V p x a(p) is the Berry curvature and B = 
V x A(x) is the 3-dimensional magnetic field. The deter- 
minant of [jab] is det[7 Q b] = (1 + fl ■ B) 2 . We see that the 
invariant phase space volume becomes \/de t [j a b ] d 3 xd 3 p , 
where ^/det [j a b] = |1 + O • B| indicates the change of 
phase space volume with time [231 ] . 

Quantum kinetic equation for Wigner functions. 
In a quantum kinetic theory, the classical phase-space 
distribution f(x,p) is replaced by the Wigner function 
W(x,p) in space-time x and 4-momentum p, defined as 
the ensemble average of the Wigner operator [25-27] for 
spin- 1/2 fermions, 



where e^ vpa is the Levi-Civita anti-symmetric tensor, 
with e 0123 = -e i23 = 1- ^(x,p) and ^(x,p) are 
the vector and axial-vector component of the Wigner 
function, which will give rise to the vector and axial- 
vector current, respectively, after integrating over four- 
momentum. 

We assume a system close to local equilibrium un- 
der a constant external field F^ v . Therefore, Y^x.p) 
and srf^ix^p) will depend on x only through fluid four- 
velocity u{x), temperature T(x), chemical potential 
and chiral chemical potential fis(x). We can determine 
the analytic form of the Wigner function in terms of 
{p, F^ v ,m,T, ^, fi 5 } from Eqs. dHEl). We further assume 
that the space-time derivative d x and the field strength 
F^ v are small variables of the same order and can be used 
as parameters in the power expansion of ^ and si^ so 
that one can use an iterative scheme to solve 1^ and ^ 
order by order. The unique forms of ^ and si^ to the 
first order that satisfy Eqs. ([6][9]) was determined in Ref. 
24j and reads, 



W, 



a/3 



d A y 
(2tt)4 
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v i/>p(x + )U(x + ,x-)il> a (x-), (5) = p»5(p 2 )Z + -p v [uV - uW 



x ± 



where ip a and ipp are Dirac spinor fields, x± 
are two space-time points centered at x with space 
time separation y, and the gauge link U(x+,x-) = 
cxp[— iQ f x+ dz p A fl (z)} ensures the gauge invariance of 

W a p. Here Q is the electromagnetic charge of the 
fermions, and A p is the electromagnetic vector poten- 
tial. Note that we use the metric convention g^ v = 
diag(l, — 1, — 1, — 1). To simplify the quantum kinetic 
equation under a background field we consider a massless 
and collisionless fermionic system in a constant external 
electromagnetic field F pv in the lab frame. Since we only 
consider an Abelian classical background field, we have 
dropped the path ordering in the gauge link. The Wigner 
function is a matrix in Dirac space and satisfies the quan- 
tum kinetic equation (H^, 7 M (p M + \V^)W{x,p) = 0, 



2"" L ~ " " 'd(u-p) 

Q PlJ [u»B» -u»B»}Z Q 5'(p 2 ) 
Qe^uxPpE^ZoS'ip 2 ), 



S(P 2 ) 



(10) 



where 2 = {"?,*(), Z Q = {V ,A ), Z = {A ,V ), with 
the first order solutions Vq and Aq given by 



[V ,A ] = 0(sU-P)[(fs,R + fs,L),(f S ,R-fs, L )] 

s=±l 

f 2 1 



(27r) 3 e s( - u -P-»x)/ T + l 



,(x = R,L), (11) 



where 7 M 's are Dirac matrices and V p = d£ — QF p v d v p 
The Wigner function should contain information about 



where R{L) denotes the right (left)-handed fermions and 
fJ-R,L = M± P>5 0- We have used notations E a = u p F api 
B a = (\l2)e ailvp u»F v P and w M = (l/2) Wff u^u CT , 
which depend on x via the fluid velocity u(x). We use 2£§ 
to denote the zero-th order term p p S(p 2 )Zo in Eq. (fT0|) 
and to denote the rest, i.e. the first order terms. 
Note that solutions of T M and gtf^ in Eq. (fTU)) are static 
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ones where (15 and T are constants, /1 = const. — QE ■ x 
and it CT (x — ut/u°) is in a solenoidal form with d ■ u = 0. 

We can derive the vector and axial-vector current 
and the energy-momentum tensor from and s^^ in 
Eq. (fTUf by integrating over momentum: j 7 ' = J d 4 p1 /l1 , 
jg = J dV^, and T» v = \ J d 4 p(p'"1 /u + p"^). The 
current jl* contains two parts proportional to magnetic 
field and vorticity, known as the CME and CVE [lfE |jj , 
respectively. So we have shown that both effects are di- 
rect consequences of the quantum kinetic equation for 
the Wigner function 24] . The axial- vector current jg also 
has two parts proportional to magnetic field and vorticity 
which can be regarded as a sort of dual effects to CME 
and CVE, respectively 0, H3|. Conservation equations 
for j M , jg (with anomaly) and can then be derived 
by taking divergences of j*, jg and [Hj]: = 0, 

d M jg = -0sE ■ B, and d„T^ = QF»>>j p . 

Relativistic chiral kinetic equation. — Now we try to 
derive a new form of relativistic chiral kinetic equation 
from Eq. ([7]). We start with Eq. j7]) for the zero-th order 
Wigner function , 



Vp3% = 



(d*-QF^d;){ P »5(p 2 )z } 

= 5(p 2 )\ P »£%-Qp»F flv d;}Z = 0, (12) 



where Z is the phase space distribution function and 
given in Eq. (jlll) . Eq. (fT"2")) is a Vlasov-like equa- 
tion, in which we can extract dx a /dr = p a /too an d 
dp^/dr = Qp^F^/mo, where r is the invariant time, 
too is a mass scale (note: it is not the fermion mass since 
we are considering massless fermions), and Qp^F^/mo 
is a general Lorentz force. It is interesting to see that 
the integration over p$ in the local rest frame gives the 
Boltzmann equation without collision terms, 



s[dt + sv • V x 

z g =±i 

+s<j(E + sv x B) • V sp ] 



(13) 



We see that the classical kinetic equation for fermions and 
anti-fermions can be derived and that the anti-fermion 
part has the same structure as the fermion one except 
the sign of charges and momenta. 

Now we consider Eq. ([7]) for the first order Wigner 
function if 3 /*, 



V^iT/ 1 = QS( P 2 )[(u ■ b)B» - (b ■ B) 
+Q5(p 2 )e^ pa u u b p E a d p Z 
~Q 2 8(p 2 )(E ■ B)b a d p J + Q8(p 2 )oj p FP°dvz 
+QS(p 2 )(p a F** - p^F^)b v uj ( dPZ , 
= (14) 

where V = —p a /p 2 . We will show that 5(p 2 )b a is a 
4-dimensional monopole in Euclidean momentum space. 



Combining Eq. (IT2l and Eq. (fl4|) . we obtain the relativis- 
tic chiral kinetic equation 



-> 6(p 2 ) 



dx a 
~~dT 





dp a 
~dr 



dl 



f R/L = 0, (15) 



where the upper/lower sign corresponds to the right/left- 
hand distribution, and dx a / dr and dp a /dr are given by 



dx a 



too— = p° ±Q[(u-b)B» -(b-B)i 



dr 



dp° 



-e^u v b p E a ], 



toq— = -Qp P FP° ±Quj p FP' 7 ±Q{p°F^ -jfiF**) 
dr 

xb v uJ S TQ 2 (E ■ B)b a . (16) 
Here we have used notations 

f R /L = l(V ±A Q ) = £ 6(su-p)f s>R/L . (17) 

Z s=±l 

In the local rest frame with u a = (1,0), B a = (0, B) and 
E a = (0, E), the Berry curvature appears in Eq. (fTB"| as 
dx°/dr = [p ± Q(h ■ B)]/too with b the spatial com- 
ponent of b a , which is similar to 1 + fl ■ B in Ref. [19j . 
Turning off the oj terms in Eq. , we obtain 



dx. 
~dr~ 
dp 
~ah 



— [p±Q(6 B + E x b)] 

to 

-^-[p E+(pxB)]±^-(E-B)b, (18) 
to too 



which is similar to Eqs. (14-15) in Ref. [19]. Part (but not 
all) of the iv terms can be generated by the replacement 
QB a — >• —ju ■ p)Lu a . The same relation is also observed 
in Ref. [19|. It is interesting to observe that there is a 
term pmQuj p F pa in in Eq. ([IB]). This shows that a 
vorticity behaves like a charged particle interacting only 
with the magnetic field but not the electric field. 

Phase space continuity equation. — We now take spa- 
tial and momentum divergence of dx a /dr and dp" /dr 
respectively as 



m d ( 



dx a 

~dV 
dp a 
~dV 



0. 



= T^(u-p)(u } -E)±^(E.B),(19) 



where we have used the conditions u p d pl ui 1 ' — u p d p B v = 
0, d„u» = 0, d p B p = 2(lu ■ E), d p u v = e purX u T uj x , 
e^UpdvEf. = 0, < 



^u v B a = 0, e 



a pa/3 



-2$ 



a/3 

[/'H 



and dlb° 



-2/p 2 . It is interesting to see that 



dx° /dr is conserved but dp a /dr is not. Using Eq. fp~5|) 
and dv[5{p 2 )] = 2b a S(p 2 ), we arrive at the phase space 
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continuity equation 

d, 



^)!r,L 



^S(p )f R/L 



±^-(u-b)(oj-E)S(p 2 ).f R/L . 
m 



(20) 



Note that there is an anomalous source term proportional 
to uj ■ E. 

Anomaly and 4~ dimensional monopole. - From 
Eq. (|10p , the anomalous conservation law of the left- and 
right-hand current can be derived as 



R/L 



T2nQ 2 (E-B) J d(u-p)S(u-p)f_ 
Q 2 



R/L 



(21) 



where j R i L — (j^ ± j£ )/2, and we have used from 
Eq. (H7J): f R/L (u-p = 0) = 1/(2tt) 3 . In deriving Eq. (JUJ, 
we completed all integrals in Minkowski space. In order 
to see the connection to the 4-dimensional monopole and 
the relativistic chiral kinetic equation, we can carry out 
the integral in Euclidean space as follows 

d»f R/L = ±Q 2 (E-B) J d 4 p6(p 2 )b°drj R/L 



±Q 2 (E ■ B)-lm 

7T 



d P— 



te p- 



-,/dllR/L 



±Q 2 (E ■ B)-lm 

7T 



d 4 p E 



iep 2 E 



Ve^TIr/l 



TQ (E ■ B)-lm 

7T 



d p E 



xd p J 



P'e 



P% 



JR/L 



(22) 



Note that the anomalous term comes from the E ■ B term 
in dp a /dr in Eq. ([T6| . We have taken analytic continu- 
ation po — ipi and p 2 = —p%. Note that in the second 
equality of Eq. (|2"2"j) the poles in Minkowski space are 
±|p| =p ie, in order to avoid these 



Po 



±y/\p\ 



ie 



poles in Wick rotation the integral limit of po should be 
[— ioo,ioo] which corresponds to [—00,00] for the J34 in- 
tegral. We have also used 

d^(p E /p E ) = 2^(p E ). (23) 

Although for p E ^ 0, we have (p E / 'p%) = 0, but the 
integral is non- vanishing, since 

J d^pEdl E {p E IVE) = j>dS^p E /p% = 2^. (24) 



Note that the n- volume of n-sphere or the hyper-surface 
area of (n+l)-ball with radius R is given by S n — 
[27r( n+1 )/ 2 /r((n + l)/2)]i? n . So we see that 8{p 2 )b a plays 
the role of a 4-dimensional monopole in Euclidean mo- 
mentum space, where the singular point of the monopole 
is located at p — 0. This is similar to the 3-dimensional 
case where the singular point is at p [3, [l9[ . 

In summary, we have shown that the Berry curva- 
ture and a 4-dimensional monopole in Euclidean momen- 
tum space emerge in a new chiral kinetic equation with 
manifest Lorentz covariance. The chiral anomaly has 
been shown to be given by the flux of this 4-dimensional 
monopole. The phase space continuity equation has also 
been derived but has an anomalous source term propor- 
tional to the contraction of vorticity and electric field. 
The derivation of the chiral kinetic equation is quite gen- 
eral and valid for relativistic fermionic system. It would 
be interesting to investigate the effect of the collision 
terms in the future. 

Note added: During the completion of this work, we 
learned that Son and Yamamoto were also working on 
the similar topic [28j]. 
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